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Abstract 

m 

We study the scaling properties of a georouting scheme in a wireless multi-hop 
network of n mobile nodes. Our aim is to increase the network capacity quasi linearly 
with n while keeping the average delay bounded. In our model, mobile nodes move 
according to an i.i.d. random walk with velocity v and transmit packets to randomly 
chosen destinations. The average packet delivery delay of our scheme is of order 1/v 
and it achieves the network capacity of order t - £ t - . This shows a practical 
throughput-delay trade-off, in particular when compared with the seminal result of 



\ Gupta and Kumar which shows network capacity of order Jn/ logn and negligible 

delay and the groundbreaking result of Grossglausser and Tse which achieves network 
capacity of order n but with an average delay of order y/n/v. We confirm the generality 
of our analytical results using simulations under various interference models. 



1 Introduction 

X 

Gupta and Kumar [lj studied the capacity of wireless networks consisting of randomly located 
nodes which are immobile. They showed that if each source node has a randomly chosen 
destination node, the useful network capacity is of order C\Jnj log n where n is the number 
of nodes and C is the nominal capacity of each node. However, if the nodes are mobile 
and follow i.i.d. ergodic motions in a square area, Grossglauser and Tse [2] showed that the 
network capacity can rise to 0(nCf\ by using the mobility of the nodes. Note that in this 
case, a source node relays its packet to a random mobile relay node which transmits this 
packet to its destination node only when they come close together, i.e., at a distance of order 
\j\Jn. Therefore, the time it takes to deliver a packet to its destination would be of order 
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1 We recall the following notation: (i) f(n) = 0(g(n)) means that there exists a constant c and an integer 
N such that f(n) < cg(n) for n > N. (ii) f(n) — Q(g(n)) means that there exists two constants c\ and C2 
and an integer N such that c\g{n) < f(n) < C2g(n) for n > N. 
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Table 1: Network Capacity vs. Delivery Delay Trade-off. 

y/nL/v where v is the average speed of the nodes and L is the length of the fixed square area 
where nodes are deployed. In contrast, in Gupta and Kumar's result [T], the packet delivery 
delay tends to be negligible, although the network capacity drops by a factor of y/n logn. 

In this article, we aim to maximize the capacity of mobile networks while keeping the mean 
packet delivery delay bounded with increasing number of nodes. For relaying packets towards 
their destinations, mobile nodes use our proposed georouting strategy, called the Constrained 
Relative Bearing (CRB) scheme. We show that, in a random walk mobility model, this 
strategy achieves a network capacity of order ; — C with a time to delivery of order 

OJ 1 J log n log log n J 

L/v. Our main contribution is summarized in Table [TJ Note that in random walk mobility 
models, nodes have free space motion and move in straight lines with constant speed. This 
mobility model is a subclass of the free space motion mobility model. Therefore, we can also 
extend our result to mobility models where the average free space distance i is non zero. 

Consider an example of an urban area network in a fixed square area of length L with 
number of nodes n = 10 6 , nominal bandwidth C = 100 kbps and delay per store-and-forward 
operation of 1 ms. The average packet delivery delay for Gupta and Kumar's case would be 
around one second but with a network capacity of 10 Mbps. In the case of Grossglauser and 
Tse, the network capacity would increase to about 100 Gbps but if the straight line crossing 
time L/v is about one hour (e.g., with cars as mobile nodes), the time to delivery would be 
around one month. However, our model of using mobility of nodes along with the proposed 
CRB scheme, to relay packets to their destinations, would lead to a network capacity of 10 
Gbps with time to delivery of about one hour. 

This article is organized as follows. We first summarize some important related works and 
results in Section [2j We discuss the models of our network and CRB scheme in Sections [3] 
and |H respectively. The analysis of capacity and delay can be found in Section [S] and we 
confirm this analysis using simulations in Section El We also discuss a few extensions of our 
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work in Section [7| and concluding remarks can be found in Section [HJ 



2 Related Works 

The main difference between the proposed models in the works of Gupta and Kumar [TJ 
and Grossglauser and Tse [2] is that in the former case, nodes are static and packets are 
transmitted between nodes like "hot potatoes", while in the latter case, nodes are mobile 
and relays are allowed to carry buffered packets while they move. Both strategies are based 
on the following model: if p n is the transmission rate of each node, i.e., the proportion of 
time each node is active and transmitting, the radius of efficient transmission is given by 
r n ~ L,/— when n approaches infinity for some constant k > which depends on the 
protocol, interference model, etc. 

In the context of [TJ, the number of relays a packet has to traverse to reach its destination 
is h n = 0(l/r n ). Consequently, np n C must be divided by h n to get the useful capacity: 
np n C/h n = 0(C ^Jp n n). In order to ensure connectivity in the network, so that every source 
is able to communicate with its randomly chosen destination, p n must satisfy the limit 
Pn < 0(l/logn). This leads to Gupta and Kumar's maximum capacity of 0(C '^/n/\ogn) 
with "hot potatoes" routing. 

In contrast, in the context of [2], the network does not need to be connected since the packets 
are mostly carried in the buffer of a mobile relay. Therefore there is no limit on p n other 
than the requirement that it must be smaller than some e < 1 that depends on the protocol 
and some other physical parameters. Thus r n is 0(1/ \/n). In Grossglauser and Tse's model, 
the source transmits the packet to the closest mobile relay or keeps it until it finds one. 
This mobile relay delivers the packet to the destination when it comes within range of the 
destination node. Such a packet delivery requires a transmission phase which also includes 
retries and acknowledgements so that the packet delivery can be eventually guaranteed. 

The proposed model of [2] requires a GPS-like positioning system and the knowledge of 
the effective range r n . The estimate of r n could be achieved via a periodic beaconing from 
every node, where each beacon contains the position coordinates of the node, so that a node 
knows the typical distance for a successful reception. However, the relay cannot rely on 
beaconing in order to detect when it is in the reception range of the destination. The reason 
is that a node stays in the reception range of another node for a short time period of order 
r n /v = lj\fn and this cannot be detected via a periodic beaconing with bounded frequency 
since p n = 0(1) (the frequency of periodic beaconing should be of 0(\/n)). We may also 
assume that the destination node is fixed and its cartesian coordinates are known by the 
mobile relay. Otherwise, if the destination node is mobile, there would be a requirement 
for this node to track its new coordinates and disseminate this information in the wireless 
network as in [31 H] . 

It is also interesting to note that Diggavi, Grossglauser, and Tse |5j showed that a constant 
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throughput per source-destination pair is feasible even with a more restricted mobility model. 
Franceschetti et al. [6] proved that there is no gap between the capacity of randomly located 
and arbitrarily located nodes. Throughput and delay trade-offs have appeared in [71 IE] 
where delay of multi-hop routing is reduced by increasing the coverage radius of each trans- 
mission, at the expense of reducing the number of simultaneous transmissions the network 
can support. We will show that, in our work, we have a delay of 0(l/v) and throughput 
per source-destination pair of 0( lognlc | glogr J . If we take the notation of a/ a(n) in [3 E] to 
measure the average distance traveled toward the destination between two consecutive emis- 
sions of the same packet, then we will show that our scheme yields \J a(n) = 0(1/ logn). If 
we compare with the result of [7J E], we should have a throughput of ®( log n J^ gn) but our 
scheme delivers a higher throughput by a factor greater than ^fn. In fact, if I is the average 
free space distance of the random walk, then our scheme yields a/ a(n) = 6 ^ i +1 1 ogn ) • The 
apparent contradiction comes from the fact that the authors in [7J [8] consider a mobility 
model based on brownian motion. This corresponds to having t = and, in this case, our 
scheme would be equivalent to the "hot potatoes" routing of [lj with a/ a(n) = Q(r n ). Let 
us point out that the brownian motion mobility is an important yet worst case model and it 
is not realistic for real world situations such as urban area mobile networks. In the section 
devoted to generalizations, we extend our result to fit a more general mobility model where 
mobile nodes follow fractal trajectories with £ = £ n = 0(1/ logn) and the throughput of our 
scheme remains of Q(-, r— ■, ). 

V log n log log n I 

On the practical side, many protocols have been proposed for wireless multi-hop networks. 
These protocols may be classified in topology-based and position-based protocols. Topology- 
based protocols [HI [lOl [UJ need to maintain information on routes potentially or currently 
in use, so they do not work effectively in environments with high frequency of topology 
changes. For this reason, there has been an increasing interest in position-based routing 
protocols. In these protocols, a node needs to know its own position, the one-hop neighbors' 
positions, and the destination node's position. These protocols do not need control packets 
to maintain link states or to update routing tables. Examples of such protocols can be found 
in [12J, [131 El E3 UHl H71 HH1 [19]. In contrast to our work, they do not analyze the trade-off 
between the capacity and the delay of the network under these protocols and their scaling 
properties. 



3 Network and Mobility Settings 

We consider a network of n mobile nodes with their initial positions uniformly distributed 
over the network area. Each mobile node transmits packets to a randomly chosen fixed node, 
called its destination node, which is also randomly located in the network area. We assume 
that mobile nodes are aware of their own cartesian coordinates, e.g., by using GPS or from 
the initial position, a mobile node could use the knowledge of its motion vector to compute 
its cartesian coordinates at any given time. 
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Initially we consider that only mobile nodes participate in the relay process to deliver packet 
to its destination node. The case where the fixed nodes may also participate in the relay 
process is discussed in Section [3 A mobile node should be aware of the cartesian coordinates 
of the destination node of a packet it carries. Indeed it can be assumed that this information 
is included in all packets or is relayed with the packets. Hence our model only requires that 
a source or relay node is aware of the cartesian coordinates of the destination node which 
is assumed fixed. Note that if the destination node is mobile, a mechanism to disseminate 
its updated cartesian coordinates in the network can be used, e.g., [U 0]. However, this is 
outside the scope of this paper as we particularly focus on the throughput-delay tradeoff. 

With the available information, a mobile relay can determine: 

- its heading vector, which is the motion vector when its speed is non zero, 

- its bearing vector, which is the vector between its position and the position of a packet's 
destination; and, 

- the relative bearing angle, i.e., the absolute angle between its heading and bearing vectors. 

In the example of Fig. [H node A is carrying a packet for node D. This figure also shows 
the heading vector of mobile node A and its bearing vector and relative bearing angle for 
destination node D. Note that a mobile relay may carry packets for multiple destinations 
but can easily determine the bearing vector and relative bearing angle for each destination 
node. 

4 Model of CRB Scheme 

In this section, we will present the parameters and specifications of the model of our georout- 
ing scheme. 

4.1 Parameters 

We define the parameters 9 C , called the carry angle, and 9 e , called the emission angle. Each 
mobile node carries a packet to its destination node as long as its relative bearing angle, 9, 
is smaller than 9 C which is strictly smaller than 7r/2. When this condition is not satisfied, 
the packet is transmitted to the next relay. 

4.2 Model Specification With Radio Range Awareness 

In the following description, we initially assume that each node is aware of the effective range 
of transmission r n . This means that there is a periodic beaconing that allows this estimate 
to be made. In Section H~3l we will investigate how to specify our model without an estimate 
of the effective range r n . 
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Mobile node 




Figure 1: Figurative representation of our model. Unfilled circles represent the potential 
mobile relays for packet transmitted by node A for node D. 

Assume that node A is carrying a packet for node D. The velocity of node A is denoted by 
v(A). 

- If node A is within range of node D, it transmits the packet to D; otherwise, 

- if the relative bearing angle is smaller than 8 C , node A continues to carry the packet; 
otherwise, 

- node A transmits the packet to a random neighbor mobile node inside the cone of angle 
6 e , with bearing vector as the axis, and then forgets the packet. 

In order to better understand the model of our georouting scheme, consider the example 
in Fig. [TJ Assume that node A is out of range of node D and, because of that, it cannot 
deliver the packet directly. Now, if 9 < 8 C , node A will continue to carry the packet for node 
D. Otherwise, it transmits the packet to one of the random mobile relays, represented by 
unfilled circles in the figure. 

4.2.1 Transmission procedure 

To transmit the packet towards another mobile node, node A shall proceed as follows: 

- it first transmits a Call-to-Receive packet containing the positions of nodes A and D; 

- a random mobile node B which receives this Call-to-Receive packet can compute the angle 
(AB,AD). If this angle is smaller than 8 e , it replies with an Accept-to-Receive packet 
containing an identifier of node B\ 

- node A sends the packet to the first mobile node which replied with an Accept-to-Receive 
packet. 

The first node which sends its Accept-to-Receive packet notifies the other receivers of the 
Call-to-Receive packet, to cancel their transmissions of Accept-to-Receive packets. There 
may be more than one (but finite) transmissions of Accept-to-Receive packets in case two or 
more receivers are at distance greater than r n from each other. 
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Note that this procedure does not need any beaconing or periodic transmission of hello 
packets. The back-off time of nodes, transmitting their Accept-to-Receive packet, can also be 
tuned in order to favor the distance or displacement towards D, depending on any additional 
optional specifications. 

4.3 Model Specification Without Radio Range Awareness 

The estimation of r n would require that the nodes employ a periodic beaconing mechanism. 
If such a mechanism is not available or desirable, the CRB scheme relies on the signal to 
interference plus noise ratio (SINR) for transmitting packets to their destinations or random 
mobile relays. In other words, a mobile node can relay a packet to its destination node or 
another mobile node only if the SINR at the receiver is above a given threshold. 

Note that in this case, the specification of the transmission procedure is also modified so 
that it terminates when the final destination receives the packet. To transmit the packet 
towards its destination node or another mobile node, node A shall proceed as follows: 

- it first transmits a Call-to-Receive packet containing the positions of nodes A and D; 

- if node D receives this packet, it responds immediately with an Accept-to-Receive packet 
with highest priority. Node A, on receiving this packet, relays the packet to node D; 
otherwise, 

- the procedure of selecting a random mobile node, as the next relay, is similar to the 
procedure described in Section 14.2.11 A random mobile node B, which receives the Call-to- 
Receive packet, computes the angle (AB, AD). If this angle is smaller than 9 e , it responds 
with an Accept-to-Receive packet; 

- node A relays the packet to the first mobile node which sent its Accept-to-Receive packet 
successfully. The first node which transmits its Accept-to-Receive packet also makes the 
other receivers to cancel their transmissions of Accept-to-Receive packets. 

5 Performance analysis 

We will show that our georouting scheme is stable as long as the average transmission rate of 
each mobile node is p n = 0(1/ log log n). We will also show that the number of transmissions 
per packet is of O(logn) and this would lead to a useful network capacity of Q(C lognl o glogn )- 

We assume that the network area is a square area and without loss of generality we assume 
that it is a square unit area. The mobile nodes move according to i.i.d. random walk: from a 
uniformly distributed initial position, the nodes move in a straight line with a certain speed 
and randomly change direction. The speed is randomly distributed in an interval [v m j n , v max ] 
with t> min > 0. To simplify the analysis, we assume that t> min = t> max = v. We also assume 
that each node changes its direction with a Poisson point process of rate r. When a mobile 
node hits the border of the network, it simply bounces like a billiard ball. This leads to the 
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isotropic property (Jacquet et al. [20]): at any given time the distribution of mobile nodes is 
uniform on the square and the speed are uniformly distributed in direction independently of 
the position in the square. 

We assume that the radius r n of efficient transmission is given by 

/log logn 

r n — \ P , 

V nn 

for some /3 > 0. Therefore, the average number of neighbors of an arbitrary node at an arbi- 
trary time is (3 log log n. In order to keep the average cumulated load finite, the nodes have 
an average transmission rate of p n = a loglogn - Therefore, the actual density of simultaneous 
transmitters is - 3 - r J2 j . 

p log log n 



5.1 Methodology 

The parameters of interest are the following: 

- The delay D n (r) of delivering a packet to the destination when the packet is generated in 
a mobile node at distance r from its destination node. 

- The average number of times F n (r) the packet changes relay before reaching its destination 
when it has been generated in a mobile node at distance r from its destination node. 

In order to exhibit the actual performance of our proposed CRB scheme, we aim to derive an 
upper-bound on the parameters D n (r) and F n (r). In the next two sub-sections, we assume 
w.l.o.g. that there is always a relay node, to receive the packet, in the emission cone (as the 
node density and angle, 9 e , are sufficiently large) when a relay change must occur. 



5.2 Delivery Delay 

In the quantity D n (r), we ignore the queueing delay which can become apparent when several 
packets could be in competition in the same relay to be transmitted at the same time. We 
analyze the delay under the hypothesis that store and forward delays are negligible (these 
delays would be negligible as long as the queue length is bounded). 

Theorem 5.1. We have the bound 

D n (r) < —^rr . (1) 

Proof. During a relay change, the new relay is closer to the destination than the previous 
relay. Ignoring relay changes that take zero time, and neglecting the distance decrement 
during relay change, the packet moves at constant speed v with a relative bearing angle 
always smaller than 6 C . ■ 



8 




Figure 2: Figurative description of relay change due to turn. At position 1, 6 < 9 C and node 
A carries the packet for node D. At position 2, node A changes its heading vector and must 
transmit the packet. 



5.3 Number of Relay Changes 

There are two events that trigger relay changes. 

1. Relay change due to turn, i. e., the mobile node, carrying the packet, changes its heading 
vector such that the relative bearing angle becomes greater than 9 C . 

2. Relay change due to pass over, i. e., the mobile node keeps its trajectory and the relative 
bearing angle becomes greater than 6 C . 

Consider a packet generated at distance r from its destination. Let F^{r) be the average 
number of relay changes due to turn. Equivalently, let F£(r) be the average number of relay 
changes due to pass over. Therefore, we have F n {r) = F^(r) + F%(r) and we expect that the 
main contribution of O(logn) in F n (r) will come from F£(r). 

5.3.1 Number of Relay Changes Due to Turn 

We prove the following theorem: 
Theorem 5.2. We have the bound 

F*(r) < r . 

6 r , V COs(8 r ) 



n \ 



Proof. We consider the case in Fig. [2] and assume that a mobile node is carrying a packet 
to its destination located at distance r. The node changes its direction with Poisson rate 
t. When the node changes its direction, it may keep a direction that stays within angle 9 C 
with the bearing vector and this will not trigger a relay change. This occurs with probability 
— . Otherwise, the packet must change relay. But the new relay may have relative bearing 
angle greater than 8 C which would result in an immediate new relay change. Therefore, at 
each direction change, there is an average of relays. Multiplied by D n (r) this gives our 
upper-bound of FHr). 
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Mobile node 
at position 




Destination node 'D' 



Figure 3: Figurative description of relay change due to pass over. At position 1, 9 < 9 C and 
node A carries the packet for node D. At position 2, node A has the same heading vector 
but 9' = 9 C and it must transmit the packet. 

Note that we have not considered the turn due to bounces on the borders of square map. 
But it is easy to see via straightforward geometric considerations that they cannot actually 
generate a relay change. ■ 

5.3.2 Number of Relay Changes Due to Pass Over 

We prove the following theorem: 
Theorem 5.3. We have the bound 

7rtan(6> c ) ( r\ 

Proof. Here we consider the case of Fig. [3] We assume that a mobile node at distance r, 
from its destination, has a relative bearing angle equal to 9. If it keeps its trajectory (i.e., 
does not turn), it will need to transmit to a new relay when it passes over the destination, 
i.e., when it arrives at a distance of p{9,r) = ^ffl , r from the destination. The function of 
9 p(9,r) is bijective from [0, 9 C ] to [0, r]. For x £ [0,r] let p~ l (x,r) be its inverse. 

Assume that r is the distance to the destination when the relay receives the packet or just 
after a turn. Thus the angle 9 is uniformly distributed on [0,# c ], i.e., with a constant 
probability density The probability density of the pass over event at x < r (assuming no 
direction change) is therefore 

1 d _ 1( . sin(# c ) tan(p -1 (x, r)) 1 

9 c dx^ ' 9 C cos(p _1 (x, r))r 9 C x 

Since p~ l (x,r) < 9 C , the point process where the packet would need a relay change due to 
pass over is upper bounded by a Poisson point process on the interval [r n , r] and of intensity 
equal to ta ~ for x £ [r n ,r]. 

Since a relay change due to pass over corresponds to an average of f- relays, and neglecting 
the decrement of distance during each transmission phase, we get 

f r 7rtan(0 c ) dx 7rtan(6U fr\ 

F " (r)= l^^- = ^^ log UJ • 
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We have thus 




7r - 6 C t 
9 C cos(9 c )v 



r + 



7rtan(# c ) 



log ( — 



Therefore we have a main contribution of O(logn) relay changes that comes from log(l/r n ). 
The result holds because we assume that there is always a receiver in each relay change. In 
the next sub-section we remove this condition to establish a result with high probability. 

5.4 Number of Relay Changes With High Probability of Success 

In the previous subsection we assumed that there is always a receiving relay in the emission 
cone at each relay change and we said that the relay change is always successful. The case 
with failed relay change would introduce additional complications. For example one could 
use the fixed relays if the packet cannot be delivered to a mobile relay. Anyhow, to simplify 
the present contribution, we will show that with high probability, i.e., with probability 
approaching one when n approaches infinity, every relay change succeeds. 

Theorem 5.4. With high probability on arbitrary packets, all relay changes succeed for this 
packet and are in average number F n {r) and the delay is D n {r). 

Proof. We use a modified stochastic system to cope with failed relay changes. The modifi- 
cation is the following: when there is no relay in the emission cone during a relay change a 
decoy mobile relay is created in the emission cone that will receive the packet. Each decoy 
relay is used only for one packet and disappear after use. Notice that the modified system is 
not a practical scheme in a practical network. The analysis in the previous section still holds 
and in particular F n (r) is now the average unconditional number of relay changes (including 
those via decoy relays) for any packet starting at distance r from destination. 

Let P n (r) be the probability that a packet starting at distance r has a failed relay change. 
The probability that a relay change fails is equal to (1 — ^ e r^) n_1 ~ e ~ n9 ^ r " = (logn)^^^. 
Therefore the average number of failed relay changes E n (r) < F n (r) (log n)~ ,/5 ''i? which tends to 
zero when (3^ > 1, since F n {r) = O(logn). The final result comes since P n (r) < E n {r). ■ 

6 Simulations 

We performed simulations with CRB georouting scheme under two contexts: 

1. a simplified context where the network is modeled under unit disk model; 

2. a realistic context where the network operates under slotted ALOHA and a realistic 
SINR interference model is considered. The simulations of CRB scheme are stressed 
to the point that the motion timings are not so large compared to slot times. 
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6.1 Under Disk Graph Model 



In this section, we consider a network of n mobile nodes. We assume that all nodes have the 
same radio range given by 



Each mobile node moves according to an i.i.d. random walk mobility model, i.e., it starts 
from a uniformly distributed initial position, moves in straight line with constant speed and 
uniformly selected direction and reflects on the borders of the square area (like billiard balls). 

In the next section (Section I6.2[) . we will further explore the effect of interference on the 
simulations, but for the moment we only consider a source mobile node and its randomly 
located destination node which is fixed. We adopt the disk graph model of interference, i.e., 
two nodes are connected or they can exchange information if the distance between them 
is smaller than a certain threshold (called radio range), otherwise, they are disconnected. 
A mobile node relays the packet only if the relative bearing angle, i.e., the absolute angle 
made by the heading vector and the bearing vector, becomes greater than 9 C . Otherwise, it 
continues to carry the packet. 

6.1.1 Simulation parameters and assumptions 

The purpose of our simulations is to verify the scaling behavior of average delay and number 
of hops per packet with increasing number of nodes in the network. Therefore, the number 
of mobile nodes, n, in the network is varied from 10000 to two million nodes. The values of 
other parameters, which remain constant, and do not impact the scaling behavior are listed 
as follows. 

(i) Parameters of CRB scheme, 9 C and 9 e , are taken to be 7r/6. 

(ii) The speed of all mobile nodes is constant, i.e., 0.005 unit distance per slot. 

(Hi) All mobile nodes change their direction according to a Poisson point process with mean 

equal to 10 slots. 
(iv) The value of constant factor (3 is assumed to be equal to 40. 

6.1.2 Results 

We have evaluated the following parameters. 

(i) Average delay per packet. 

(ii) Average number of hops per packet. 

We considered the Monte Carlo Method with 100 simulations. The delay of a packet is 
computed from the time when its processing started at its source mobile node until it reaches 
its destination node. Figure H] shows the average delay per packet with an increasing number 
of nodes. We notice that as n increases, the average delay per packet appears to approach a 
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Figure 5: Average number of hops per packet. 



constant upper bound which can be computed from (ITl) . Figure shows the average number 
of hops per packet with increasing values of n. 



6.2 With slotted ALOHA under SINR interference model 

In this section, we will present the simulations of CRB georouting scheme with a transmission 
model which does not rely on the estimate of r n and is based on the required minimal SINR 
threshold. 



6.2.1 Transmission model 



Our transmission model is as follows. Let Pi be the transmit power of node i and 7^ be the 
channel gain from node i to node j such that the received power at node j is Pffij- The 
transmission from node i to node j is successful only if the following condition is satisfied 

P ™ > K , 



where K is the desired minimum SINR threshold for successfully receiving the packet at the 
destination and Nq is the background noise power. For now, we ignore multi-path fading or 
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shadowing effects and assume that the channel gain from node i to node j is given by 

1 



I 



a 



where a > 2 is the attenuation coefficient and z% is the location of node i. 
6.2.2 Simulations under SINR interference model 

For the theoretical analysis in Section [51 we have assumed that the effective range of successful 
transmission is 

^loglogn 

which requires that the mobile nodes have an average transmission rate of p n — (3/ log log n. 
In other words, if the mobile nodes emit packets at the given average rate, the average 
distance of successful transmission under SINR interference model is of 0(r n ) and the results 
from theoretical analysis are applicable as well. 

We assume that time is slotted and mobile nodes determine their relative bearing angles at 
the beginning of a slot. We also assume that all nodes are synchronized and simultaneous 
transmitters in each slot emit a Call-to-Receive packet at the beginning of the slot. Moreover, 
we also assume that fixed nodes do not emit any packet except, maybe, an Accept-to-Receive 
packet in response to a transmission by a mobile node. 

In our simulation environment, n mobile nodes start from a uniformly distributed initial 
position and move independently in straight lines and in randomly selected directions. They 
also change their direction randomly at a rate which is a Poisson point process. Each mobile 
node sends packets towards a unique destination (fixed) node, and all destinations nodes are 
also uniformly distributed in the network area. 

In order to keep load in the network finite, the packet generation rate at a node, p n , should 
be of 0(p n /X n ) where X n is the expected number of transmissions per packet. From the 
theoretical analysis, we know that 

X n = o(\og (^X) +c, 



where c is a constant if 8 C is non-varying. In our simulations under SINR interference model, 
we assume that the knowledge of r n is not available and mobile nodes use minimal SINR 
threshold for successfully receiving a packet. We also assume that each mobile node generate 
packets, destined for its unique fixed destination node, at a uniform rate given by 

1 

Pn ~ /3 1 log(^)loglogn ' (2) 

for some fi\ > and (3 2 > 0. We ignored the value of constant c and have observed that the 
simulation results are asymptotically correct because, with n increasing, value of c should 
be insignificant as compared to the 0(log(n/ fa)) factor. 
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6.2.3 Simulation parameters and assumptions 

The purpose of our simulations is to verify the scaling properties of network capacity, delay 
and number of transmissions per packet with increasing number of nodes in the network. 
The number of mobile nodes, n, in the network is varied from 250 nodes to 100,000 nodes. 
All nodes use uniform unit nominal transmit power and the background noise power N Q is 
assumed to be negligible. The values of other parameters are listed as follows. 

(i) Parameters of CRB scheme, 8 C and 9 e , are taken to be 7r/6. 

(ii) The speed of all mobile nodes is constant, i.e., 0.01 unit distance per slot. 

(Hi) All mobile nodes change their direction independently and randomly according to a 
Poisson point process with mean equal to 10 slots. 

(iv) The values of constant factors (5\ and f3 2 are assumed equal to 500 and 1 respectively. 

(v) SINR threshold, K, is assumed equal to 1. 

(vi) Attenuation coefficient, a, is assumed equal to 2.5. 

In our simulations, we make the following assumptions. 

(i) Each mobile node generates an infinite number of packets, at rate p n , for its respective 
destination node. 

(ii) A mobile node may carry, in its buffer, its own packets as well as the packets relayed 
from other mobile nodes. Therefore, it may have more than one packet in its buffer which 
it must transmit because their respective relative bearing angles are greater than 9 C . In 
such a case, it first transmits the packet which is furthest from its destination. 



6.2.4 Results 

We have examined the following parameters. 

(i) Throughput capacity per node, A n . 

(ii) Average number of hops, h n , and transmission attempts, t n , per packet. 
(Hi) Average delay per packet. 

The throughput capacity per node, A n , is the average number of packets arriving at their 
destinations per slot per mobile node. With n increasing, throughput capacity per node 
should follow the following relation 



n ftlog^loglogn ' 

for some < rj < 1 which depends on K, a and protocol parameters. Note that the values 
of these constants do not affect the asymptotic behavior of A n which is also observed in our 
simulation results. 
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Figure 6: Verification of network throughput capacity with plots of m\ and m p . 




Figure 7: Simulated (solid lines) and theoretical (dotted lines) network throughput capacity, 
nX n , and network packet generation rate, np n . 



In order to verify the asymptotic character of simulated packet generation rate and through- 
put capacity, we have analyzed the parameters m p and m\ which are given by 

m p = p n ( fk log ( -=- J log log(n 



0: 



m x = A n ( fa log ( — ) loglog(n) 



0: 



From the definition of p n in <^j, the value of m p should be constant at 1 whereas, with n 
increasing, value of m\ should converge to the constant r\. From Fig. [6j value of r] is found 
to be approximately equal to 0.45. Figure [7] shows the simulated and theoretical packet 
generation rate, np n , and throughput capacity, n\ n , in the network. The theoretical values 
of np n and nX n are computed from (j5J) and (jSJ). 

Figure M shows the average number of hops, h n , and transmission attempts, t n , per packet. 
The value of t n is slightly higher than the value of h n because of the possibility that a suc- 
cessful receiver may not be found in each transmission phase, i.e., in the cone of transmission 
formed with 9 e . With n increasing, h n and t n are expected to grow in 0(log(n//?2))- To 
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Figure 8: Average number of hops, h n , and transmission attempts, t n , per packet. 




Figure 9: Verification of number of hops and transmission attempts with plots of and 
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Figure 10: Average delay per packet. 



verify this character in simulation results, we examine the parameters m^ and m t given by 

m h = K- — — — , 
log(£) 
1 

m t = t r , 



If the values of h n and t n are in 0(log(n//?2)), the values of and mt should approach a 
constant value which is the case in Fig. |9j 
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The delay of a packet is computed from the time when its processing started at its source 
mobile node until the time it arrives at its destination node. Figure [10] shows the average 
delay per packet. As the number of mobile nodes increase, the average delay appears to 
approach a constant value. 

It can be observed that when n is small, the average number of hops per packet is almost 
of 0(1) which also means that the average delay per packet is of 0(1) and the network 
throughput capacity is of 0(rjn): although, in simulation results, it is bounded by the network 
packet generation rate which is of O( logral " glogra ). This can be observed in Fig. [TJ [51 and [Till 
The reason is that when n is small, the number of simultaneous transmissions in the network 
is also small and packets can be delivered by the mobile nodes, directly to their destination 
nodes, in 0(1) hops. As n increases, number of simultaneous transmitters increase and 
consequently the effective transmission range of each transmitter shrinks. Therefore, the 
dominant factor in the number of transmissions per packet comes from the fact that a mobile 
relay has to be close to the destination, to deliver a packet. According to theoretical analysis, 
h n and t n grow in O(logn) which is also observed in the simulation results. Simulations also 
show that, asymptotically, network throughput capacity is of 0( lognl " glogw ) and average 
delay per packet is of 0(l/t>) which complies with our theoretical analysis. 



7 Extensions and general mobility models 

In our discussion, we primarily focussed on the capacity-delay tradeoff and thus for the initial 
sake of clarity assumed that the fixed nodes can only receive packets destined for them. We 
could also consider a slight variation in the specification of the model of CRB scheme such 
that the fixed nodes also participate in the routing of packets to their destination nodes. For 
example, during a transmission phase, if a packet cannot be transmitted to its destination 
node or relayed to a random mobile neighbor in the cone of transmission, it can be relayed to 
a fixed node. This fixed node must emit this packet immediately to its destination node or to 
any mobile relay in the neighborhood. Note that this will also help increase the connectivity 
of the network. 

The condition about i.i.d. random walks can be relaxed and the result about the expected 
number of relay changes will still be valid. In other words, the i.i.d. random walk model 
can be seen as a worst case compared to realistic mobility models. If the mobile relays move 
like cars in an urban area, then we can expect that their mobility model will significantly 
depart from the random walk. Indeed cars move toward physical destinations and in their 
journey on the streets toward their destination, their heading after each turn is positively 
correlated with the heading before the turn. This implies that the probability that a relay 
change is needed after a turn is smaller than it would be under a random walk model, where 
headings before and after turn are not correlated. Furthermore on a street, the headings are 
positively correlated (consider Manhattan one-way streets) and in this case a relay change 
due to a pass over will have more chances to arrive on a relay with good heading (one half 
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Figure 11: Illustration of self-similar trajectories in urban areas. 

instead of 9 c /n). Again this would lead to less relay changes due to pass over. 

The result still holds if we assume that the turn rate r depends on n and r = r n = O(logn). 
In this case, the mobility model would fit even better for the realistic mobility of an urban 
area. Indeed the trajectories of cars should be fractal or self-similar, showing more frequent 
turns when cars are close to their physical destination (different than packet destination) or 
when leaving their parking lot. In this case, the overall turn rate tends to be in O(logn) 
with a coefficient depending on the Hurst parameter of the trajectory. This would lead to 
the same estimate of O(logra) relay change per packet. 

Figure [11] illustrates a self-similar trajectory in an urban area. It shows a two-dimensional 
trajectory (upper half) and its traveled distance (lower half). The successive turns are 
indicated by T 1; . . . ,T 7 . The trajectory after any turn Tj looks like a reduced copy of the 
original trajectory. The CRB scheme may need some adaptation to cope with some unusual 
street configurations, e.g., to replace the cartesian distance with the Manhattan distance in 
the street map. 



8 Conclusions 

We have examined asymptotic capacity and delay in mobile networks with a georouting 
scheme, called CRB, for communication between source and destination nodes. Our results 
show that CRB allows to achieve the network capacity of 0(-, — ) with packet delivery 

r J v log n log log n ' L J 

delay of 0(1) and transmissions per packet of O(logn). It is noticeable that this scheme 
does not need any sophisticated overhead for implementation. However, in this case, the 
mobile nodes must be aware of their position via a GPS system, for example. 

We have shown the asymptotic performance via analytical analysis under a unit disk graph 
model with random i.i.d. walks. The analytical results have been confirmed by simulations 
and in particular under ALOHA with SINR interference model. We have seen that the 
performance of CRB can be maintained even with non i.i.d. random walks, the latter being 
a worst case scenario. However, this latter result would require that the mobile nodes stay 
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within same heading for 0(1/ logn) time. A next step would be to analyze the performance 
of this scheme on real traffic traces in urban areas. 
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